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Abstract 

Background: Human mesenchymal stromal cells from the bone marrow (BMSCs) are widely used as experimental 
regenerative treatment of ischemic heart disease, and the first clinical trials using adipose-derived stromal cells (ASCs) 
are currently being conducted. Regenerative mechanisms of BMSCs and ASCs are manifold and in vitro pretreatment 
of the cells with growth factors has been applied to potentially enhance these properties. When characterizing the 
transcriptional activity of these cellular mechanisms in vitro it is important to consider the effect of the growth factor 
treatment on reference genes (RGs) for the normalization of qPCR data. 

Results: BMSCs and ASCs were stimulated with vascular endothelial growth factor A- 165 (VEGF) for one week, and 
compared with un-stimulated cells from the same donor. The stability of nine RGs through VEGF treatment as 
well as the donor variation was assessed using the GenEx software with the subprograms geNorm and Normfinder. 
The procedure of stepwise elimination was validated by poor performance of eliminated RGs in a normalization 
experiment using vWF as target gene. Normfinder found the TATA box binding protein (TBP) to be the most stable 
single RG for both BMSCs and ASCs. The optimal number of RGs for ASCs was two, and the lowest variance for vWF 
normalization was found using TBP and YWHAZ. For BMSCs, the optimal number of RGs was four, while the two-RG 
combination producing the most similar results was TBP and YWHAZ. 

Conclusions: A common reference gene, TBP, was found to be the most stable standalone gene, while TBP and 
YWHAZ were found to be the best two-RG combination for qPCR analyses for both BMSCs and ASCs through the 
VEGF stimulation. The presented stepwise elimination procedure was validated, while we found the final 
normalization experiment to be essential. 

Keywords: Adipose-derived stromal cell, ASC, ADSC, MIQE, qPCR, Reference gene, Mesenchymal stromal cell, MSC, 
Vascular endothelial growth factor, VEGF 



Background 

Human bone marrow derived stromal cells (BMSCs) and 
adipose-derived stromal cells (ASCs) have made a sig- 
nificant impact on the field of regenerative medicine. 
Within the field of cardiac regeneration, BMSCs have 
been used frequently in clinical trials while the first clin- 
ical trials with ASCs stimulated with vascular endothelial 
growth factor (VEGF) is currently being conducted 
by our group [1,2]. BMSCs are the most extensively 
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investigated cell type, however adipose tissue has re- 
ceived increasing interest as a source of stem cells due to 
easier harvest, a larger yield of stem cells and the fact that 
ASCs possess higher proliferative capacity compared to 
BMSCs [3]. The rationale behind the potential benefit of 
using the two cell types has primarily been suggested on 
the basis of their immunoregulatory properties, paracrine 
profiles, and capability to differentiate into various cell 
types [4-9] . When investigating these properties and how 
to enhance them in vitro, the first change within the cells 
is found at the transcriptional level. This change is de- 
tectable by quantitative real time polymerase chain reac- 
tion (qPCR) analyses, which is the preferred method for 
detecting changes in gene expression [10]. However, the 
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qPCR method is sensitive to several errors in the experi- 
mental setup that may affect results. These errors could 
occur at different stages throughout the experiment and 
encompass variations in the amount of starting material, 
nuclear extraction, RNA integrity, cDNA loading, reverse 
transcription, and qPCR efficiency [11]. Thus, in order to 
compare qPCR results from different runs, or samples 
obtained from different tissues at various time points, a 
normalization step is necessary [12]. 

The most widely used method for normalizing qPCR 
data is comparison to values from reference genes (RGs) 
analyzed simultaneously on the same material [13]. The 
abundance of the normalization candidate gene must be 
closely related to the total amount of mRNA in the sam- 
ple, in order to be able to correct for experimental errors 
[14]. In addition, it must exhibit stable expression during 
experimental conditions to be useful for normalization. 
Often an invariant endogenous gene is used as a RG 
with a constant expression level across all samples and 
experimental treatments. If the RGs are not stable 
through the treatment it could obscure actual changes 
and produce artifacts that contribute to misleading 
results and incorrect conclusions, thereby leading to un- 
reliable publications [13]. The attention regarding selec- 
tion of suitable RGs has increased in the past decade, 
from the geometric mean normalization proposed by 
Vandesompele et al, to experiment-specific RG investi- 
gations, and ultimately the "Minimum Information for 
Publication of Quantitative Real-Time PCR Experiments" 
(MIQE) guidelines [14,15]. The guidelines function as a 
universal checklist that investigators can follow, to produce 
the most reliable data when using qPCR [10]. 

In the field of cardiac regeneration, it has been pro- 
posed that endothelium-committed cells greatly enhance 
cell therapy outcome [16]. Ongoing clinical trials per- 
formed by our group serve to differentiate ASCs and 
BMSCs toward endothelial lineage for enhancement 
of angiogenic potential by stimulation with vascular 
endothelial growth factor (VEGF) [2,17,18]. This pro- 
angiogenic growth factor is known to promote endothe- 
lial cell proliferation, migration, and vessel permeability, 
and to differentiate BMSCs toward endothelial lineage 
[19,20]. The impact of VEGF treatment on ASCs and 
BMSCs is tested in vitro by the use of qPCR, where it is 
very important for both cell types to have established 
RGs that are stable through the VEGF treatment. A suit- 
able RG in ASCs is not necessarily stable in BMSCs, 
therefore, it is important to validate whether the chosen 
RGs can be used in both cell sources under specific 
experimental conditions. No studies comparing choice of 
RGs between ASCs and BMSCs have been conducted 
previously to the authors' knowledge, which makes this 
study particularly interesting with regard to the ongoing 
debate of differences and similarities between these two 



cell types. Through this study we provide a simple step- 
wise selection procedure for choosing the RGs with focus 
on the impact on the final outcome. For the purpose we 
use specially designed software GenEx, more specifically 
the algorithms Normfinder and geNorm [14,21]. 

Results 

To be able to apply the qPCR on VEGF pre-conditioned 
human mesenchymal stromal cells from the bone mar- 
row and adipose tissue, a panel of RGs (Table 1) from 
relevant literature was chosen to find the correct RG for 
normalization. After confirming RNA purity, RNA integ- 
rity, primer efficiency, and products (Figure 1), a stepwise 
elimination of candidate RGs was performed in order to 
select the best-suited RG(s) from a list of acceptable ones. 

RNA and primer quality 

RNA purity was validated by excluding protein contamin- 
ation with absorbance ratios of 1.8-2.2 at 260 nm/280 nm. 
RNA integrity was confirmed by RIN values of 10 for all 
donors (Additional file 1: Table SI). Primers (Table 1) for 
the qPCR were validated by dissociation curves (Additional 

file 1: Figure SI) and by efficiencies (E = ^10 5/o ^ e -l^ TOO) 

of 100% ±10% with a correlation coefficient R 2 between 
0.985 and 1.0 (Additional file 1: Table S2). 

RG candidate selection 

Step 7: RG candidate expression levels 

The RG candidate Cq values were generally between 15 
and 25 cycles, with no great difference between the 
VEGF treatment group and the control group for either 
ASCs or BMSCs (Figure 2). The selected RGs should 
preferably be expressed in amounts close to those of the 
target genes. The raw data for the outliers provided us 
with Cq values that were substantially lower than that of 
the gene of interest, as depicted in Figure 2 presenting 
median Cq values. Reference genes with a substantially 
greater difference in abundance compared to the gene of 
interest are not suitable. It would be incorrect to per- 
form an analysis on putative reference genes that are 
unsuitable to be considered. Here our cutoff Cq value 
difference is 10 cycles. In order not to reduce the initial 
panel of reference genes completely, we didn't exclude 
more than the three genes mentioned at this stage. 

We have investigated the impact of VEGF treatment 
in terms of initiation of endothelial differentiation in 
ASCs and BMSCs, after which endothelial target genes 
such as vWF are expressed in relatively low amounts, 
and therefore detected around 30 cycles. Since EFl-a, 
ACTB, and 18S rRNA were the most abundant RGs in 
both ASCs and BMSCs, with Cq values far from our tar- 
get, they were removed from the panel (Figure 2). This 
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Table 1 Reference genes for qPCR 



Gene Gene 



NCBI reference 
sequence 



Forward sequence 



Cellular function 



18SrRNA 18S ribosomal RNA 



ACTB beta-actin 



GUSB beta-glucuronidase 



PPIA 



RPL13 ribosomal protein LI 3a 



TBP TATA-binding protein 



NR_003286.2 



NM_001 101.3 



EF1a elongation factor- 1 alpha NM_001402.5 



GAPDH glyceraldehyde 3-phosphate NM_002046.4 
dehydrogenase 



NM_000181.3 



peptidyl prolyl isomerase A NM_021 130.3 



NM_01 2423.3 



NM_003 194.4 



YWHAZ tyrosine 3-monooxygenase/ NM_001 135702.1 
tryptophan 5-monooxygenase 
activation protein, zeta 
polypeptide 



F = 5'- 
R = 5'- 



R = 5'- 
F = 5'- 
R = 5'- 

F = 5- 
R = 5'- 
F = 5'- 
R = 5'- 



R = 5'- 
F = 5'- 
R = 5'- 



-GTAACCCGTOAACCCCA^-3' 

-CCATCCAATCGGTAGTAGCG-3' 

-CC^mGTCCCCCAAOTGA-3' 

-TGGCTGCCTCCACCCA-3' 

■AGGTGATOTCCTGAACCATCC-3' 

-AAAGGTGGATAGTCTGAGAAGC-3' 

-CAACGGAmGGTCGTATOG-3' 

•GCAACAATATCCACmACCAGAGTOA-3' 

-CTCAmGGAA^GCCGA^-3' 

-CCGAGTGAAGATCCCC^mA-3' 

-TCCTGGCATOTGTCCATG-3' 

-CCATCCAACCACTCAGTC1TG-3' 

-CATAGGAAGCTGGGAGCAAG-3' 

-GCCCTCCAATCAGTC1TCTG-3' 

-TGCACAGGAGCCAAGAGTGAA-3' 

-CACATCACAGCTCCCCACCA-3' 

-AC^GGTACATOTGGOTCAA-3' 

-CCGCCAGGACAAACCAGTAT-3' 



Small subunit of cytoplasmic 
ribosomes 

Cytoskeleton structural protein; 
motility, cytokinesis 

Translation 



Oxido-reductase in glycolysis and 
gluconeogenesis, transcription 
activation, apoptosis 

Catalyzing hydrolysis of 
B-D-glucuronic acid 

Protein folding 



Structural component of 60S 
ribosomal subunit 

RNA polymerase II transcription 
factor 

Signal transduction 



Information for the panel of selected reference genes with full name, 
cellular function. 



NCBI reference sequence, forward and reverse primer sequences, and short description of 



initial exclusion resulted in the six RGs expressed in 
closest abundance to the target gene remaining for 
further analyses. 

Step 2: Initial standard deviation-values from Normfinder 

Normfinder selects the most suitable RG based on the 
standard deviations (SDs) for intragroup donor variation 
and intergroup treatment variation. This analysis was 
performed taking the group classification of treatment 
into account, in order to verify that the candidate RGs 



were stably expressed and not regulated by treatment. 
The RG candidates with a large bias, an intergroup SD 
reaching or exceeding 0.2 were excluded (Table 2), leav- 
ing TBP, PPIA, YWHAZ, and GUSB for ASCs and TBP, 
PPIA, YWHAZ, GUSB, and GAPDH for BMSCs for the 
final analysis. 

Step 3: Final output from geNorm and Normfinder 

Due to the selection in step 2, the software should not 
regard the effect of treatment in the next step of the 




150 65 234 71 80 87 156 131 93 



Figure 1 Performance of primer amplification. qPCR endpoint products were separated by agarose gel electrophoresis. Target genes are 
indicated above the amplicon. 
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Figure 2 Expression level dot plots of RG candidates. Distribution of median Cq values for ASCs (A) and BMSCs (B) from eight donors in 
control medium (blue) and VEGF-treatment medium (red). ASCs; adipose-derived stromal cells, BMSCs; bone marrow-derived stromal cells, RG; 
reference gene, VEGF; vascular endothelial growth factor. N = 8 except for vWF where N = 3. 



selection. This means that the final GenEx comparisons 
encompass both intragroup donor variation and inter- 
group treatment variation. According to Normfinder 
ranking order, TBP was found to be the most stable RG 
candidate through the VEGF treatment for both ASCs 
and BMSCs (Figure 3A and 3B). The optimal number of 
RGs was two for ASCs and four for BMSCs based on 
the accumulated SDs, which were reduced slightly by 
0.009 and 0.003 respectively when compared to the use 
of a single RG (Figure 3C and 3D). In order to determine 
which RGs to use in combination we utilize intergroup 
treatment variation data calculated by Normfinder using 
only the remaining genes from Step 3 (Figure 3E). The 

Table 2 Step two of elimination, showing ASC and BMSC 
reference gene intergroup variations 



Gene name ASC (SD) BMSC (SD) 



GAPDH 


±0.3855* 


±0.0780* 


GUSB 


±0.0048* 


±0.0832* 


PPIA 


±0.0977* 


±0.1028* 


RPL13 


±0.2843* 


±0.2865* 


TBP 


±0.0090* 


±0.0016* 


YWHAZ 


±0.1832* 


±0.0435* 



Difference in mean expression of each gene between treated and untreated 
groups is shown as standard deviations (SD) in unit cycles. Genes were excluded 
at a cutoff value above 0.2. Those with a SD below the threshold value (*) were 
included considered for further processing. 18S rRNA: 18S ribosomal RNA, ACTB: 
beta-actin, EF1-a: Elongation factor 1 -alpha, GAPDH: gluco phosphate 
dehydrogenase, GUSB: Beta-glucuronidase, PPIA: Peptidyl prolyl isomerase A, 
RPL13a: Ribosomal protein L13-alpha, TBP: TATA box binding protein, YWHAZ: 
Tyrosine 3/tryptophan 5-monooxygenase activation protein. Genes not indicated 
with a star (*) were deactivated in GenEx. 



sum of the different possible combinations were calcu- 
lated and compared. The chosen combinations of TBP 
and YWHAZ for ASCs and TBP, GAPDH, PPIA and 
GUSB for BMSCs yielded a RG variation sum closest to 
zero. 

GeNorm uses pair-wise comparison and sequential 
elimination to select the pair of RGs with highest correl- 
ation, indicated by the average expression stability (M). 
TBP and GUSB were found to be the best correlated RGs 
for ASCs and PPIA and GAPDH for BMSCs (Figure 4A 
and 4B) with a cut-off level at 0.5. 

RG verification 

We compared the effect of normalization to various RG 
choices on the expression of a gene of interest (GOI), 
vWF, to validate our procedure of stepwise elimination 
of RG candidates and to attain a final selection of the 
most suitable RGs for future VEGF stimulation experi- 
ments. vWF was chosen as positive control since it has 
been shown to be upregulated in other experiments by 
the same treatment, namely serum-deprivation com- 
bined with VEGF stimulation [19,22]. vWT was normal- 
ized to RGs from each selection step of our procedure. 
For both ASCs and BMSCs, EFl-a was chosen from step 1 
while RPLI3 was chosen from step 2. From step 3, the sin- 
gle best RG found and the optimal number of RGs found 
by Normfinder and the two best correlated RGs given by 
geNorm were evaluated. Geometric mean of RG combina- 
tions were used, as suggested by Vandesompele et al [14]. 
For ASCs, vWF was normalized to TBP, the geometric 
mean of TBP and YWHAZ (Normfinder), or the geometric 
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Figure 3 Determination of most stable gene and optimal number of RGs. In the third step of RG selection, the Normfinder algorithm was 
used on the genes remaining following the first and second eliminations. Gene ranking according to stability through VEGF treatment of ASCs 
(A) and BMSCs (B). Accumulated SDs indicate the best number of RGs for ASCs (C) and BMSCs (D). The intergroup treatment variation as 
calculated by Normfinder during Step 3 of the RG selection process (E). The shaded histograms present the single most stable reference gene. 
ASCs; adipose-derived stromal cells, BMSC; bone marrow-derived stromal cells, SD; standard deviation. N = 8. 



ASC 



B 



BMSC 



TBP 



GUSB 



PPIA GAPDH 



Figure 4 Determination of the best correlated genes. Average expression stability measure (M) calculated by pairwise comparison of RGs 
with sequential elimination of least stable RG by geNorm for ASCs (A) and BMSCs (B). The shaded histograms indicate the two most correlated 
genes. ASCs; adipose-derived stromal cells, BMSC; bone marrow-derived stromal cells, SD; standard deviation. N = 8. 
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mean of TBP and GUSB (geNorm) (Figure 5). For BMSCs, 
the normalization was also performed to TBP alone, to 
the geometric mean of TBP, PPIA, GAPDH, and GUSB 
(Normfinder), to the geometric mean of PPIA and GAPDH 
(geNorm) (Figure 6A). Additionally, vWF was normalized 
to the geometric mean of TBP and YWHAZ on the basis 
of Figure 6B. 

Poor results with normalization to eliminated RGs 

For both ASCs and BMSCs, EFl-a was chosen from step 
1 while RPL13 was chosen from step 2 and the single 
best RG found by Normfinder, TBP, was chosen from 
step 3. When comparing the fold difference in vWF ex- 
pression for ASCs, no significant increase was observed 
for VEGF-treated ASCs compared to controls, regardless 
of the RGs used for normalization. The trend in a slight 
upregulation of vWFwas not investigated in depth, since 
this was not the aim of the study. For assessing our elim- 
ination procedure, comparison of the use of single RGs 
gave different results for the previously eliminated RGs 
compared to TBP, The variance in RPL13 results was 
markedly larger than the variance in TBP results, while 
EFl-a results showed a tendency of down-regulation of 
vWF, in contrast to results from all other normalizations. 
For BMSC, normalizing to any of the RG combinations 
or TBP alone included in step 3 of selection, proved 
vWF to be significantly increased by the VEGF treat- 
ment. The results are not significant when normalizing 



to EFla or RPL13, and their variance was greater than 
that for TBP. 

Evaluating combinations of RGs 

Comparing the results from step 3 normalizations for 
ASCs, TBP alone produced results almost identical to ei- 
ther of the two combinations. Variance was lowest in 
TBP and highest in the geNorm combination, yet still 
almost identical (Figure 5). 

For BMSCs, the fold expression increase obtained 
using TBP normalization was slightly below that for 
the Normfinder combination of RGs, while the geNorm 
combination resulted in a higher fold expression level 
(Figure 6A). Since the geNorm combination is the two 
best correlated RGs and the Normfinder combination 
was the best calculated number of RGs, we trust the ex- 
pression level found using the Normfinder combination 
to be the most precise. 

We further investigated whether a combination of 
two RGs could replace the four RG combination sug- 
gested by Normfinder for practical purposes. We found 
the geometric mean of TBP and YWHAZ to be the 
combination with the closest related and nearly identi- 
cal vWF fold change level, also with the lowest SEM 
across samples (Figure 6B). This combination was com- 
pared to the others in Figure 6A. When comparing the 
results of vWF expression between the different RG 
combinations, only TBP is significantly different from 
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Figure 5 Difference in vWF expression normalized to different RGs. Normalization of vWF to various RGs for ASCs. The RGs were chosen 
from geNorm and Normfinder results. Two previously eliminated RGs, EF1-a and RPL13 are included. ASCs; adipose-derived stromal cells, RG; 
reference gene, vWF: von Willebrand Factor. N = 3 and error bars represent SEM. 
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Figure 6 Comparison of BMSC vWF expression results from two-RG-combinations. Normalization of vWF to various RGs for BMSCs (A). The 
RGs were chosen by the same method as those in Figure 5. "denotes significant (p < 0.05) vWF change, and t denotes significant correlation with 
results using the Normfinder combination of four-RGs. Test of two-RG combinations relative to results from the Normfinder combination of four 
genes (B). 2" AACq data were calculated for each RG combination using the geometric mean of the selected RGs, and compared to the mean of 
the results using Normfinder in the RG verification step. § denotes the combination producing results most closely related to those obtained using 
the Normfinder four-RG combination, and the combination producing the lowest standard deviation. BMSCs; bone marrow-derived stromal cells, 
vWF; von Willebran Factor, RG; Reference Gene, C; Combination. N = 3 and error bars represent SEM. 



the Normfinder combinations. However, only the results 
using the combination of TBP and YWHAZ, found from 
C4 in Figure 6B, were significantly correlated with the 
Normfinder combination. 

Discussion 

We compared different RGs for qPCR to detect the most 
optimal RGs for human ASCs and BMSCs, as these cell 
types are currently used in clinical trials, where they are 
stimulated with VEGF before injection into the myocar- 
dium of patients with chronic ischemic heart disease. 
When investigating the effect of the growth factor on 
the cells on transcriptional level, it is necessary to per- 
form RG optimization for the production of properly 
normalized, reliable, and reproducible results. 



The best single RG for both cell types in our experi- 
mental setup was TBP, which has been found to be 
stable through several other interventions using different 
stem cell lines or cancer cells [15,23,24]. TBP and 
YWHAZ were found to be the best two-RG combination 
for qPCR analyses of VEGF stimulation with serum 
deprivation in both ASCs and BMSCs. 

Only four RGs were included in the final analysis for 
ASCs. It is noteworthy that the reduced panel from 
which an RG or combination of RGs was to be selected, 
were the exact same RGs that Fink et al have shown to 
have stable expression in ASCs treated with hypoxia 
[15]. The other candidate RGs were excluded on the 
basis of having intergroup variation above the 0.2 cycle 
cutoff or because their expression level was significantly 
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higher in comparison to potential genes of interest. This 
excluded RPL13 which was found by several groups to 
be the most stable RG in mesenchymal stem cells 
through different treatments [25,26]. 

Though Normfinder found the optimal number of 
RGs to be two (Figure 3C), and geNorm found GUSB to 
be the RG best correlated with TBP (Figure 4A), the dif- 
ference between using TBP alone or in combination with 
GUSB was minimal. The inclusion of GUSB only leads 
to a reduction of 0.009 in accumulated SD (Figure 3C), 
and the variation in fold change of vWF between the 
two normalizations was insignificant (Figure 5). Simi- 
larly, the combination of TBP and YWHAZ, calculated 
from Normfinder, would give the same results as TBP 
alone. The geometric mean of this combination even 
had a slightly lower variance compared with that of TBP 
and GUSB. We suggest that TBP can be used alone or in 
combination with YWHAZ or GUSB for ASCs in this 
experimental setup. 

For BMSCs, five RGs were included in the final ana- 
lysis. The elimination explained earlier lead to the exclu- 
sion of EFl-a, which is considered to be a very stable 
RG, and the inclusion of GAPDH, which has been shown 
to be rather unstable through different treatments 
[15,27,28]. GeNorm identified the best combination of 
two genes to be PPIA and GAPDH (Figure 4B). However, 
the use of these two RGs leads to a higher accumulated 
SD (Figure 3D), and their geometric mean produced a 
larger fold increase in vWF expression compared to the 
other combinations (Figure 6A). 

The fact that normalizing to the geometric mean of 
the two genes chosen by geNorm produced a larger fold 
increase in vWF compared to the other combinations is 
an example of geNorm calculating pair-wise correlation 
rather than stability when selecting the best suited RGs. 
The number of RGs to be used as suggested by Norm- 
finder was the combination of four genes. The use of 
two RGs would increase the accumulated SD by 0.026 
(Figure 3D), but due to practical reasons we further ana- 
lyzed combinations of two genes to find the normalization 
that would produce the most similar result for vWF 
expression (Figure 6B). Since TBP and YWHAZ produced 
results closest to the combination of the four RGs calcu- 
lated by Normfinder and displayed the lowest SEM, we 
obtain a combination of RGs similar to that for ASCs, 
which could be useful in experimental setups. Except for 
the use of TBP as single RG, none of the vWF expression 
results using the various combinations were significantly 
different from those using the four-RG Normfinder com- 
bination. However, the use of the other combinations 
could still produce Type 1 or Type 2 errors due to skewed 
expression levels or high degree of variance. The C4 com- 
bination from Figure 6B, TBP and YWHAZ, was the 
only combination significantly correlated to the four-RG 



Normfinder combination. This result is of course depen- 
dent on our choice of the four-RG Normfinder combin- 
ation as the most stable. In contrast to BMSCs, the 
different RG combinations for ASCs produced more 
similar vWF expression results. This could be due to the 
seemingly larger impact of the VEGF treatment on 
BMSCs, which could likewise have larger influence on 
the RG expression, compared to standard culture. We 
have recently investigated the effect of VEGF treatment 
with serum-deprivation on ASCs, and found that the 
treatment did not produce endothelial differentiation to 
the same extent as shown earlier on BMSCs stimulated 
similarly [19,22]. 

The importance of using a final GOI-specific test as 
supplement to results from the software before choosing 
RGs is especially underscored by the BMSC results. 
Simply following the results from geNorm would have 
resulted in an overestimation of the VEGF effect, which 
could severely bias an entire study. The software is very 
useful in the elimination process towards the GOI- 
specific test, which is evident by our results with RGs 
eliminated in the earlier steps. The RG candidate EFl-a 
was eliminated early on in the selection process with an 
expression level far from that of vWF. For ASCs, using 
EFl-a as an RG would be misleading while in BMSCs 
there is a large variance. Normalizing to RPL13 also re- 
sults in a large variance for both cell types as expected, 
since it was eliminated early on for high intergroup vari- 
ation. The fact that EFl-a and RPL13 would be unsuit- 
able to use as RGs supports our RG selection procedure. 

Conclusions 

For reliable and reproducible qPCR analyses, selection of 
the correct RGs for a given experimental setup is of 
prime importance. We have presented a step-by-step 
procedure for selecting RGs, with emphasis on a follow- 
up test with a gene of interest as the final step in the 
selection process. Genes were eliminated based on high 
intergroup variations, after which Normfinder and geNorm 
were applied for further selection, with a final comparison 
of vWF expression levels for optimal RG selection. We 
found that the combination of TBP and YWHAZ as RGs is 
stable for normalization of VEGF-stimulated ASCs and 
BMSCs. 

Methods 

Experimental design 

ASCs were isolated from lipoaspirate obtained from 
eight healthy donors of which one was male and seven 
were female (age between 28-57 years; mean age 43.4). 
BMSCs were obtained from bone marrow aspirate from 
eight healthy donors, two males and six females (age 
20-46 years; mean age 27.8). The use of ASCs and 
BMSCs from healthy volunteers was approved by the 
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National Ethical Committee protocol no. H-3-2009-119. 
All donors agreed to and signed the informed consent. 
After isolation and cultivation, ASCs and BMSCs were 
both stimulated with VEGF and compared with untreated 
control. 

Bone marrow preparation and BMSC isolation 

50 ml bone marrow aspirate was obtained from the iliac 
crest by needle aspiration under local anesthesia. The 
sample was diluted 1:2 with phosphate-buffered saline 
pH 7.4 (PBS, -Ca 2+ -Mg 2+ , Gibco, Invitrogen, Denmark, 
cat.no. 10010-015). Mononuclear cells (MNCs) were 
harvested by gradient centrifugation on Lymphoprep 
(1077 g/cm 3 , Medinor, Denmark, cat.no. 1114547), 
washed with PBS and counted using NucleoCounter® 
NC-100™ (Chemometec, Denmark) according to man- 
ufacturers instructions. Primary cell cultures of MNCs 
were established by seeding 2 x 10 7 cells/T75-flask (Nunc, 
Thermo Scientific, Denmark, cat.no. 156494) in complete 
medium containing Dulbeccos Modified Eagle Medium, 
low glucose (1 g/1) (DMEM) supplemented with 25 mM 
HEPES and L-Glutamin, (PAA Laboratories, Austria, cat. 
no. E15-808), 10% Fetal Bovine Serum Farma grade (FBS, 
PAA Laboratories, cat.no. Al 1-512) and 1% Penicillin/ 
Streptomycin (Gibco, cat.no. 15140-122). The cells were 
incubated in standard conditions at 37°C in humid air with 
5% C0 2 . The medium was changed 5 days after initial 
seeding, and subsequently every 3-4 days. 

Lipoaspirate preparation and ASC isolation 

Approximately 100 ml lipoaspirate was obtained from 
liposuctions of subcutaneous abdominal fat performed 
under local anesthesia. The lipoaspirate was washed twice 
with PBS to remove residual blood. The adipose tissue 
was digested by incubation with collagenase (Collagenase 
NB4 (0.6 PZ U/ml, Serva GmbH, Germany) dissolved 
in HBSS (+CaCl 2 + MgCl 2 , GIBCO, cat.no. 14065-049) 
diluted to a concentration of 2 mM Ca 2+ ) at 37°C for 
45 min. under constant rotation. The collagenase was neu- 
tralized with complete medium and the suspension was 
filtered through a 100 \A mesh (Cell Strainer, BD Falcon, 
cat.no. 352360); centrifuged at 1200 g for 10 min. at room 
temperature, re-suspended and counted. MNCs were 
seeded at a density of 4.5xl0 6 cells/T75-flask in complete 
medium and incubated at standard conditions. After two 
days in culture, cells were washed with PBS to remove 
non-adhering leukocytes, and complete medium added 
anew with subsequent change of media every three-four 
days. 

Cell culture 

When the culture reached a confluence level of approxi- 
mately 90%, cells were washed with PBS, detached with 
3 ml TrypLE® (TrypLE® Select, Gibco, cat.no. 12563-029) 



for 10 min. at 37°C and neutralized with 7 ml complete 
medium. The suspension was centrifuged at 300 g for 
5 min. at room temperature, counted, and seeded at 
3.5xl0 5 cells/T75-flask for the experimental setup. The 
cells exhibited stem cell characteristics, by adhering to 
plastic, expressing stem cell markers and being able to 
differentiate [22,29-31]. 

Stimulation with VEGF 

BMSCs and ASCs were cultured in complete medium 
until their confluence was estimated to be 80%. They 
were then either kept in complete medium (control), or 
changed to serum-deprived medium (DMEM + 2%FBS + 
196P/S) added 50 ng/ml recombinant human VEGF-A 165 
(rhVEGFA 165 , R&D Systems, USA, cat.no. 293-VE-CF). 
All media was renewed every two-three days and cells 
were cultured for one week after which they were har- 
vested for further processing. 

Nucleic acid extraction 

Total RNA was extracted using Qiagen RNeasy® Mini 
Kit (Qiagen Hamburg GmbH, Hamburg, Germany, 
cat.no. 74106) according to the manufactures protocol. 
1 ml syringe (Omnific-F 1 ml, B.Braun Melsunger AG, 
Germany, cat.no. 300013) was used to homogenize the 
lysed cells before applying the Qiagen protocol. RNA pur- 
ity was measured using NanoDrop® 1000 Spectrophotom- 
eter (Thermo Scientific, MA, USA), and the eluate was 
stored at -80°C. RNA purity and integrity were confirmed 
using RNA Nano Chips (Agilent, Santa Clara, CA, cat.no. 
5067-1521) and the Agilent 2100 Bioanalyzer by following 
instructions of the Agilent RNA 6000 Nano Kit. 

Reverse transcription 

The cDNA synthesis reaction was prepared using Affinity- 
Script (Stratagene, Denmark, cat.no. 600559) in an eight- 
tube strip (0.1 mL, Micro Amp™, Appied Biosystems®, 
Invitrogen, cat.no. 4358293) on ice. The total reaction 
volume was 20 [A with 0.5 [ig RNA, 10 [A cDNA synthesis 
master mix, 3 \i\ Oligo (dT) primer, 1 \A AffinityScript 
RT RNase block enzyme mixture, and RNAse-DNAse 
free water to 20 \A total volume. Tubes were closed with 
caps (Applied Biosystems, cat.no. N801-0535) and the 
reactions were performed with an initial stage of 25°C for 
5 min., 42°C for 45 min. and 95°C for 5 min. (Veriti 96 
well fast thermal cycler, Applied Biosystems model no. 
9901). Following synthesis, the cDNA was stored in 
aliquots at -20°C. 

Quantiative real-time PCR 

qPCR was performed in triplicate per donor for each 
group within the same qPCR-run. A calibration curve 
was run simultaneously on the RG candidate tested on 
the donors. Only data obtained from runs fulfilling the 
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same criteria for efficiency and correlation coefficient as 
the primer verification was included for analysis. The Cq 
threshold was set to the value of 0.1 for all qPCR runs. 
Brilliant II SYBR'green QPCR Low ROX master mix 
(Agilent, cat.no. 600806) was used with a total reaction 
volume of 25 ul in 96-well optical reaction plates (Agilent, 
cat.no. 401333) with 5 ul of diluted cDNA. The plate was 
sealed with optical plastic caps (Agilent, cat.no. 401425). 
qPCR was performed using Mx3000 (Stratagene) and the 
results were collected using Mx3000 version 4.0 software 
for Windows (Stratagene). The reaction was initiated by 
heating to 95°C for 10 min., followed by 40 cycles elong- 
ation at 60°C for 1 min. and denaturation at 95°C for 
30 sec. 

To verify the chosen RGs, a normalization experiment 
was set up. vWF was used as target gene, and normal- 
ized to different combinations of RG candidates. The 
level and the standard deviations of the fold changes 
between VEGF treated cells and controls were compared 
for normalization to different RG combinations. The 
fold changes in vWF expression between VEGF treat- 
ment and controls were calculated with the 2" AACq - 
method. 

Gel electrophoresis 

A 3% 3-1 NuSieve agarose gel (Lonza, Switzerland, cat. 
no. 50090) was made according to the laboratory proto- 
col with lx TAE buffer diluted from 50x TAE buffer 
(Qiagen, cat.no. 129237). For visualization of qPCR 
products, 10% non-toxic GelStar Nucleic Acid Gel Stain 
(Lonza, cat.no. 50535) was added to the melted agarose 
before pouring the solution into a plastic well with 
combs. After gel solidification the combs were removed 
and the solid gel transferred to an electrophoresis tank 
(BioRad, CA, USA) and covered with 1 x TAE buffer. For 
every 8 ul PCR product, 2 ul Gelpilot 5x loading dye 
(Qiagen, cat.no. 1037650) was added and the 10 ul sample 
is loaded by directly injecting into the wells. A negative 
control was included to show no contamination of the 
product. Gelpilot 1 kb ladder (Qiagen, cat.no. SM0318) 
were used. The gel was run between 70 and 100 volt in 
45 min.- 60 min. and visualized under ultraviolet light. 

Data analysis 

Step 1. Cq values were plotted in Microsoft Excel and 
assessed visually for outliers in the lower end. 

For further analysis, we used the GenEx software 
(MultiD Analyses AB, Sweden), and more specifically its 
algorithms Normfinder and geNorm. In short, Normfin- 
der selects the most suitable RG based on the standard 
deviations (SDs) for intragroup and intergroup variation, 
while geNorm compares the RGs pair-wise and performs 
sequential elimination resulting in the best correlated 



RGs. Average Cq-values of the technical triplicates were 
used for analysis. 

Step 2. For each cell type, the mean Cq-values of all 
RG candidates from all donors were entered in GenEx 
for the treatment and control groups. 

Step 3. The function Normfinder was performed 
taking group classification of treatment into account, in 
order to verify that the candidate RGs were stably 
expressed and not regulated by treatment. The RG 
candidates with SDs reaching or exceeding 0.2 were 
removed from the panel by inactivation in the Data 
Manager of GenEx before proceeding to the next step. 
Normfinder was run again, without taking groups into 
account to give the best estimate of the genes' stabilities 
and an assessment of the optimal number of RGs. Lastly, 
geNorm was run to determine the best correlated RGs. 
The AACq-method was used when assessing the effect 
of VEGF treatment on the expression of vWF normal- 
ized to various RGs, with subsequent use of paired t-test 
and Pearsons Correlation in IBM SPSS. Results were 
considered significant at ^-values below 0.05. 

Additional file 



Additional file 1: Table SI. RNA Quality. Table S2. Primer Quality. 

Figure SI. Dissociation Curves. 
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